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Sensitivity of atomic transition frequencies to variation of the fine structure constant a — /he 
increases proportional to where Z is the nuclear charge. Recently several lines of heavy ions 
Ge II, Sn II and Pb II have been detected in quasar absorption spectra. We have performed accurate 
many-body calculations of a^-dependence of transition frequencies (g-coefficients) for these atoms 
and found an order of magnitude increase in sensitivity in comparison with atomic transitions 
which were previously used to search for temporal and spatial variation of a in quasar absorption 
systems. An interesting feature in Pb II is highly non-linear dependence on o? due to the level 
pseudo-crossings. 

PACS numbers: PACS: 31.25.-v, 31.25.Eb, 31.25.Jf 



I. INTRODUCTION 

Theories unifying gravity with other interactions sug- 
gest a possibility of temporal and spatial variation of fun- 
damental constants of Nature. A review of these theories 
and measurement results can be found in Q]. A very 
sensitive many-multiplet (MM) method to search for the 
variation of the fine structure constant a = /he by 
comparison of quasar absorption spectra with laboratory 
spectra has been suggested in ref 0|- It requires calcu- 
lations of relativistic corrections to atomic energy levels 
which allow one to find dependence of atomic transition 
frequencies on a: 

uo = bJa + qx (1) 

where x = c? /oi^ — 1 « 25a /a. Here and ao are the 
laboratory values, u and a are the rest value of transition 
frequency and the fine structure constant for an atom or 
ion in a remote cloud of ionized gas located on a distance 
up to 12 billion light years from us. Calculations of q for 
a number of atomic transitions have been performed in 
Refs. [10, H 0. Webb et al. 0, S H |l3 used MM 
method [3 and found possible evidence of a variation, 
while other groups 11, 12] have used the same method 
but found no indications of the variation ( Note, how- 
ever, that works 0,11,001 are based on data from Keck 
tele scop e located in Northern hemisphere while works 
[HI US used data from Very Large Telescope located in 
Southern hemisphere. Therefore, the disagreement be- 
tween the results, in principle, may be attributed to the 
spatial variation of a.). 

To improve these results it would be very important to 
include into data analysis new atomic transitions where 
possible effects of a variation are larger. Indeed, now 
systematic effects are at most comparable to the ob- 
served frequency shifts corresponding to value of 5a/ a — 
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(-0.54 ± 0.12) • IQ-^ reported in 0. If transitions with 
an order of magnitude larger values of q will be used these 
systematic effects will be negligible in comparison with 
the frequency shifts produced by 5a/a ~ —0.5 • 10^'^. 

Coefficients q are small in light atoms and rapidly in- 
crease (~ Z'^) with nuclear charge Z. Recently tran- 
sitions in heavy ions Ge II {Z ~ 32, rest wavelength 
1602.4863 A), Sn II [Z = 50, rest wavelength 1400.450 
A) and Pb II {Z = 82, rest wavelength 1433.9056 A) have 
been observed in quasar absorption spectra 0. There- 
fore, in the present work we perform calculations of the 
dependence of transition frequencies on a in Ge II, Sn II 
and Pb II. Coefficients q in Pb II are an order of magni- 
tude larger than in lighter elements {Z < 31) which were 
studied in Refs. E II E IH El ■ Potentially, this 
may give an order of magnitude increase of sensitivity to 
a variation and help to eliminate possible disagreement 
between the results of different groups. 

Note that only El transitions from the ground states 
are observed in quasar absorption spectra. 

II. CALCULATIONS 

Calculations for Go II, Sn II and Pb II are difficult due 
to strong configuration mixing and level pseudo-crossing. 
As it was pointed out in our previous works energies of 
different atomic states considered as functions of may 
have very different slopes and come to crossing at partic- 
ular values of a 3. Since there can be no real crossing 
of levels with the same conserving quantum numbers the 
phenomenon is called level pseudo-crossing. Configura- 
tions are strongly mixed and the slope of curve repre- 
senting E{a'^) change very quickly in the vicinity of level 
pseudo-crossing. If this happens near physical value of a 
{a « ao) the results for the coefficients q (see Eq. (^) 
become unstable since small error in a position of level 
crossing results in a large error in q. One way of dealing 
with this problem is by fitting the experimental values 
of the Lande g-factors When states of different to- 
tal angular momentum L and total spin S are mixed 
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the resulting values of ^-factors depend on configuration 
mixing similar to g-coefficients {q li^ii 9 ~ J^di^i 
where Wi is the weight of a state i). Therefore, correcting 
the mixing coefiicients to fit the experimental g-factors 
helps to find true values of the q-coefficients as well. How- 
ever, in the case of ions considered in present work, there 
is strong configuration mixing between states of the same 
L, S and J (J is the total momentum). These states are 
the ^Dj states of the ns^nd and nsnp^ configurations 
(n = 4 for Ge II, n = 5 for Sn II and n = 6 for Pb II). The 
mixed states have the same values of (/-factors and result- 
ing g-factor does not depend on configuration mixing. In 
contrast, the values of q-coefficients are still sensitive to 
configuration mixing. 

In principle, it is possible to use experimental hyperfine 
structure (hfs) to fit the configuration mixing coefficients. 
Since single-electron matrix elements of the hyperfine in- 
teraction have very different values for s, p and d states, 
different configurations must have different values of the 
many-electron hfs matrix element. Therefore, fitting the 
experimental hfs would have the same effect as fitting of 
experimental (7-factors. Unfortunately, no experimental 
data on hfs is available for Ge II, Sn II and Pb II. 

This leaves us with energies being the only control of 
the accuracy of calculation of the q-coefficients (with ex- 
ception of few levels of Pb II where ^-factors are also 
useful). Calculations need to be done to very high ac- 
curacy for the energies to be a reliable control. The cri- 
terion is that deviation of the calculated energies from 
the experimental values must be much smaller than the 
experimental energy interval between mixed states. The 
interval between different states of Ge II, Sn II and 
Pb II is a little larger than 10000 cm~^. Therefore, the 
deviation of the calculated energies from the experiment 
should be less than ^ 1000 cm~^ or < 1% of the excita- 
tion energy from the ground state. 

A method which can produce the results of desirable 
accuracy was suggested in Ref. |0. Calculations are 
done in the V^~^ approximation. This means that the 
self-consistent Hartree Fock procedure is done for the 
quadruply charged positive ion. As it has been demon- 
strated in Ref. |1^ , removal of s and p valence electrons 
doesn't really affect the atomic core (apart from singe- 
electron energies) and V^-^ approximation is a good 
approximation for all ions with number of valence elec- 
trons ranges from 1 to 5. These include neutral atom and 
negative ion. The main advantage of this approximation 
is that core- valence correlations can be relatively easy in- 
cluded beyond the second-order of the perturbation the- 
ory. As it has been demonstrated in Ref 14] inclusion 
of higher-order core- valence correlations can significantly 
improve the accuracy of calculations. 

The ions of our current interest have three valence elec- 
trons and the effective Hamiltonian for the three-electron 
wave function of valence electrons has the form 

3 3 
^Y^hu+Y^h^j, (2) 

2 = 1 i^j 



hi{ri) is the one-electron part of the Hamiltonian 

hi = cap + (/? - l)mc2 - — + V^-^ + El . (3) 

r 

El is the single-electron operator (correlation potential) 
which describes correlations between a particular valence 
electron and core electrons. It is the same operator which 
is used for atoms and ions with one external electron 
above closed shells (see, e.g. 0,0,0)- /12 is the two- 
electron part of the Hamiltonian 

h2^-. r-FE2(ri,r2), (4) 

|ri -rsl 

E2 is the two-electron part of core-valence correlations. 
It represents screening of Coulomb interaction between 
valence electrons by core electrons. 

Note that removing E from the effective Hamiltonian 
reduces it to the effective Hamiltonian of the standard 
configuration interaction (CI) method. Since we use 
many-body perturbation theory to calculate E the tech- 
nique we use can be called the CI-I-MBPT method [l8l |. 

The details of calculations for positive ions of Ge, Sn 
and Pb will be presented elsewhere. Below we discuss 
specifics of calculations for Ge II, Sn II and Pb II. 



A. Ge II 

The Ge II ion is the lightest of three ions {Z=32) 
and the easiest from computational point of view. The 
core- valence correlations are relatively small due to small 
number of electrons in the core. The relativistic cor- 
rections are small too. The latter means in particular 
that fine structure is small and energy multiplets stay 
well apart from each other and there is no level pseudo- 
crossing. Also, there is almost no mixing between states 
of different total angular momentum L and/or different 
total spin S. The ^153/2,5/2 states of the As-^id and 4s4p^ 
configurations are still strongly mixed. However, due to 
high accuracy of the calculations the final results are very 
stable. 

We calculate Ei and E2 for the effective Hamilto- 
nian ((2Jl in the second order of the MBPT. Inclusion 
of El brings single-electron energies of Ge IV to agree- 
ment with the experiment on the level of 0.1%. No 
higher-order core-valence correlations need to be in- 
cluded. The final results are presented in Table To 
calculate g-coefficients (see Eq. (Q) we perform calcu- 
lations of the energies for x = —0.1,0,-1-0.1. Than 
q+ = 10(u;(0.1) - w(0)),g_ = 10(w(0) - w(-O.l)) and 
finally q — {q+ + q-)/2. We need both q+ and q^ to 
check for non-linear behavior of energies which is usually 
an indicator of the level pseudo crossing. For Ge II (7+ 
and g_ are practically identical for all levels considered. 
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TABLE L Energy levels and relativistic energy shift of Ge II 
(cm-i). 





State 


Enerf 


jies 


1 






Expt." 


Calculations 




-'^1/2 


0.0 










2 po 

-f^3/2 


1767.1 


1797 


-1863 




'Pl/2 


51575.5 


51512 


-3098 




^^3/2 


52290.5 


52241 


-3901 




-P5/2 


53366.7 


53342 


-5001 


4s^5s 


'S'l/2 


62402.4 


62870 


-664 


4s4p^ 


-D3/2 


65015.0 


65313 


-3387 




'^5/2 


65184.1 


65494 


-3612 


4s^5p 


2 po 


79006.2 


79386 


-940 




2 po 
-f^3/2 


79365.8 


79750 


-1317 


4s^4d 


-D3/2 


80836.1 


81444 


-2236 




'A/2 


81011.8 


81625 


-2466 
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TABLE II: Energy levels and relativistic energy shift of Sn II 
(cm-i). 



State Energies q 







Expt." 


Calculations 




2 ryo 

M/2 


0.0 










2 po 
M/2 


4251.4 


4222 


-4680 




*A/2 


46464.2 


46661 


-5930 




*P3/2 


48368.0 


48556 


-8343 




*P5/2 


50730.0 


50915 


-10558 




'S'l/2 


56885.9 


56707 


-2679 


5s5p^ 


^^3/2 


58843.8 


58806 


-7361 




^^5/2 


59463.4 


59419 


-8325 


5s^5d 


D3/2 


71405.6 


71140 


-5904 




'1)5/2 


72047.6 


71804 


-7015 


5s^6p 


2 po 
M/2 


71494.3 


71182 


-3119 




2 po 
-f^3/2 


72377.3 


72061 


-4071 



"Moore, [l^ 



B. Sn II 

The Sn II ion {Z = 50) is very similar to the Ge II 
ion. However, correlations and relativistic corrections 
are larger. It has some implication on the calculation 
scheme. It turns out that inclusion of the higher-order 
core- valence correlations does lead to significant improve- 
ment of the results for both Sn IV and Sn II ions. We in- 
clude screening of Coulomb interaction and hole-particle 
interaction in all orders of the MBPT. It is done exactly 
the same way as in our calculations for single-valence- 
electron atoms (see, e.g. Qjj). The S2 operator is still 
calculated in the second order of the MBPT. 

Final results are presented in Table |n] The q- 
coefficients are calculated exactly the same way as for 
Ge II. The resulting accuracy for energies is very good 
and the values of q-s are very stable. 



C. Pb II 

The Pb II ion {Z = 82) is the most difficult for cal- 
culations. Correlations are strong and relativistic effects 
are large too. Strong L — S interaction lead to intersec- 
tion of the fine-structure multiplets. Also, states of the 
same total momentum J are strongly mixed regardless 
of the values of L and S assigned to them. The breaking 
of the L — S scheme can be easily seen e.g. by compar- 
ing experimental values of the Lande g-factors with the 
non-relativistic values. 

We have done one more step for Pb II to further im- 
prove the accuracy of calculations as compared to the 
scheme used for Sn II. We have introduced the scaling 
factors before Si to fit the energies of Pb IV. These en- 
ergies are found by solving Hartree-Fock-like equations 
for the states of external electron of Pb IV in the v'^~* 
potential of the atomic core 

{Ho + e„)^n = 0. (5) 

Here Hq is the Hartree Fock Hamiltonian. Si is the all- 
order correlation potential operator similar to what is 
used for Sn II. Inclusion of Si takes into account the 
effect of the core-valence correlations on both the en- 
ergies (e„) and the wave functions (V'n) of the valence 
states producing the so-called Brueckner orbitals. The 
difference between Brueckner and experimental energies 
of the 4s, 4p and Ad states of Pb IV are on the level 
of 0.2 - 0.4%. To further improve the energies we re- 
place Si by /Si with rescaling factor / chosen to fit the 
energies exactly. Then the same rescaled operator /Si 
is used for the Pb II ion. It turns out that only small 
rescaling is needed. Maximum deviation of the rescal- 
ing factor from unity is 10%: /(4s) = 0.935, /(4pi/2) = 
1.084, /(4p3/2) = 1.1, /(4d3/2) = 1.07, /(4d5/2) = 1.07. 

The rescaling of the Si operator has similar effect on 
both ions Pb IV and Pb II (and Pb HI) improving agree- 
ment with experiment in all cases. Although, the effect 
on the energies of Pb II is small, it is important due to 
the level crossing and strong configuration mixing which 
can make the results to be very unstable if the accuracy 
of calculations is not high enough. 

Note that we fit the energies of Pb IV but not Pb II. 
The calculations for Pb II can be still considered as pure 
ab initio calculations since no experimental information 
about Pb II is used. 

The results are presented in Table ITTll We also present 
on Fig. I the behavior of Pb II energy levels as functions 
of (a/ao)^. Dotted lines correspond to states with J = 
1/2, short dash lines correspond to J — 3/2 and long 
dash lines correspond to J = 5/2 (the longer the dash the 
larger the J). Horizontal dashes on the right represent 
experimental energies. 

The resulting picture is very complicated. There are 
at least three level pseudo-crossing. One is for the 
and ''P5/2 states at about (a/ao)^ = 0.8, another is for 
the '*Pi/2 and ^6*1/2 states for (a/ao)^ > 1 and the third 
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TABLE III; Energy levels and relativistic energy shift of Pb II 
(cm-i). 



State Energies q 







Expt." 


Calculations 




z r}o 

^1/2 













2 no 
^3/2 


14081 


13897 


-19210 




'Pl/2 


57911 


58052 


-21000 




'P3/2 


66124 


66222 


-34500 




IPb/2 


73905 


73749 


-40700 




Sl/2 


59448 


59204 


-6800 




lDs,2 


68964 


69256 


-43900 




D3/2 


69740 


69002 


-14600 




2 po 

-^1/2 


74459 


73878 


-8180 




2 po 
^3/2 


77272 


76666 


-12240 




D3/2 


83083 


83196 


-45200 




'A5/2 


88972 


88800 


-46700 



"Moore, ^ 



one is for the and states at (a/ao)^ > 1. In 

fact, due to the level pseudo-crossing, these states are so 
strongly mixed that their assignment to particular L ~ S 
multiplets and configurations is ambiguous. 

To improve the accuracy in q we used experimental 
information about the g-factors and energy intervals be- 
tween the strongly mixed levels. These values are re- 
produced with high accuracy in the interval of (a/ao)^ 
between 0.9 and 1. Also, this interval is far from the 
positions of the level pseudo-crossings which happen at 
(a/ao) > 1 and (a/ao)^ ~ 0.8. Therefore, the results are 
stable. In this situation we use g_ as the best estimate 
of the g- values: q = q- = 10((jj(0) — w(— 0.1)) {oj = oj{x), 
where x = (a/ao)^ ~ !)■ A deviation of g = (g_ + q+)/2 
from q_ gives us a reasonable accuracy estimate in this 
case, ~ 10%. Note again that the energies are calculated 
with much higher accuracy. 

III. CONCLUSION 

We used combination of the many-body perturbation 
theory and configuration interaction method to calcu- 
late dependence of the energy levels on the fine-structure 
constant (q-coefficients) . This dependence is due to the 
Dirac relativistic corrections (which are strongly modified 
by the many-body effects). As we found in our previous 
works the Breit interaction and QED radiative correc- 
tions are not important for the atoms of interest. 

Calculations in Ge II and Sn II are straightforward 
and dependence of the transition frequencies on is 
close to the linear one {q+ « g_). The accuracy of the 
energy level calculation in Ge II and Sn II is better than 
1%, for the fine structure intervals the accuracy is better 
than 3%. Therefore, the errors in q should not exceed 



3%. The levels where q are anomalously small may be 
an exception. For example, a conservative estimate of 
accuracy for the level with ui = 62402 cm~^ in Ge II 
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FIG. 1: Energies of lowest even-parity multiplets of Pb II as 
functions of c? . 



where q — —664 cm^^ is better than 10%. The value of 
q for this level (g = —607 cm~^) was also calculated in 
our previous work Pj. The present calculation is more 
accurate. 

The case of Pb II is more complicated due to the level 
pseudo-crossing (as functions of a). The dependence on 
c? is highly non-linear (see Fig 1). However, after taking 
into account the experimental information about the g- 
factors and intervals between energy levels the errors in 
g- values are reduced to ~ 10%. 

The values of q in Sn II and especially Pb II are much 
larger than in the elements which were previously used 
to search for a variation where q ~ 1000 cm^^. For the 
detected transitions Ge II (Z = 32, rest wavelength 
1602.4863 A) q= -660 cm"!, Sn II {Z = 50, rest wave- 
length 1400.450 A) g = -5900 cm"!, and Pb II {Z = 82, 
rest wavelength 1433.9056 A) g = -14600 cm'^. The 
largest (in absolute value) calculated g-coefficient in Pb 
II is (7 « —45000 cm^^. This may give significant increase 
in sensitivity to the variation of a. 
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